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ABSTRACT 

Context. Stochasticity and physical parameter degeneracy problems complicate the derivation of the parameters (age, mass, and 
extinction) of unresolved star clusters when using broad-band photometry. 

Aims. We develop a method to simulate stochasticity and degeneracies, and to investigate their influence on the accuracy of derived 
physical parameters. Then we apply it to star cluster samples of M31 and M33 galaxies. 

Methods. Age, mass and extinction of observed star clusters are derived by comparing their broad-band UBVRI integrated magnitudes 
to the magnitudes of a large grid of star cluster models with fixed metallicity Z = 0.008. Masses of stars for a cluster model are 
randomly sampled from the initial mass function. Models of star clusters from the model grid, which have all of their magnitudes 
located within 3 observational errors from the magnitudes of the observed cluster, are selected for the computation of its age, mass, 
and extinction. 

Results. In the case of the M3 1 galaxy, the extinction range is wide and the age-extinction degeneracy is strong for a fraction of its 
clusters. Because of a narrower extinction range, the age-extinction degeneracy is weaker for the M33 clusters. By using artificial 
cluster sample, we show that age-extinction degeneracy can be reduced significantly if the range of intrinsic extinction within the host 
galaxy is narrow. 

Key words, galaxies: star clusters: general 
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1. Introduction 

Star clusters are important objects for understanding the forma- 
tion and evolution of their host galaxie s. It is consid ered that 
most of star formation is clustered ( |Lada & Lad a 2003 ), there- 
fore knowledge of the physical parameters of a cluster popu- 
lation (e.g., age, mass, chemical composition, and extinction) is 
essential for constraining the star formation history of the galaxy. 

The commonly used method for deriving the physical param- 
eters of unresolved star clusters is based on comparing of their 
integrated broad-band photometry colors to the colors of simple 
stellar population (SSP) models; see, e.g., An ders et al.[ (2004) 
and |Bridzius et al.| ([2008). However, this method is strongly bi- 
ased by the presence of two major problems: 

- degeneracy between various physical parameters of star clus- 
ters, see, e.g., |Worthey| ( |1994) , |Bridzius et al.| ( |2008| ). For ex- 
ample, a young cluster possessing high extinction can have 
colors similar to an older object without extinction - an age- 
extinction degeneracy; 

- stochasticity, which is due to the random presence of a few 
bright stars, which dominate the integrated photomet r y of 
unresolved clusters; see, e.g.,|Santos & Frogel|([l997|), |De- 



days preferred approach (e.g. ,|Popescu & Hanson|( 2009, 2 010]); 
Fouesneau & Lancon (2010 ); Fouesneau et al. (20 12j ); |Asa , D &| 
Hanson ( 20121)) is to use a Monte-Carlo sampling of stellar IMF 
to model integrated colors of clusters and build an extensive grid 
of models, to cover all possible age, mass and extinction ranges. 
Physical parameters of star clusters are then derived by compar- 
ing observations to the grid of models. 



Recently |Asa'D & Hanson| ( |2012| ) have derived the age 



veikis et al.| ( [2008] ), |Maiz Apellaniz|p009| ).^Young clusters 



can have supergiant stars that significantly redden their inte- 
grated colors. By using the SSP method, a much older age 
would be determined for these clusters. 



Although |Cervino & Luridiana] ( |2006| ) have attempted to 
describe the problem of stochasticity analytically, the nowa- 



and extinction of star clusters in the Large Magellanic Cloud 
(LMC) using the method of |Popescu & Hanson| ( |2010| ) that 
takes stochasticity into account, and compared results to previ- 
ous derivations based on an isochrone fit to the resolved color- 
magnitude diagrams (CMDs). They managed to constrain the 
age of clusters similar to the values found by the isochrone fit 
only when previously known extinctions of individual clusters 
were used. 

In this paper, a method of deriving physical parameters is 
developed and applied to star cluster samples (integrated multi- 
band photometry) in two Local Group galaxies: 1) M31, using 
a catalog by Vanseviciu s et al.| ( |2009| ), who derived cluster pa- 
rameters using SS P meth od, and 2) M33, usin g objects common 
to catalogs of Ma ( 2012 ) for photometry and San Roman et al. 
(2009), who observed in resolved conditions with the Hubble 
Space Telescope (HST) and estimated age, mass, and extinction 
based on an isochrone fit to cluster CMDs. 

Using this data and artificial cluster simulations, we demon- 
strate that when the intrinsic range of extinction within the host 
galaxy is rather narrow, it is not necessary to know the exact 
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value of the extinction for individual clusters to derive their 
physical parameters reliably. 

The paper is organized into the following sections. Section[2] 
introduces our method of deriving physical parameters of clus- 
ters when stochasticity is taken into account, §[3] describes a grid 
of simulated cluster models, § HI presents tests of the method on 



artificial cluster samples, and 
and M33 star clusters. 



applies the method to the M3 1 



2. Method of deriving age, mass, and extinction of 
star cluster 

There are presently two main methods used to derive physical 
parameters of star clusters based on a 3D grid (age, mass and 
extinction) of models, which take stochasticity effects into ac- 
count. The first is a fast^ 2 minimization approach used by, e.g., 
|Popescu & Hanson] ( |2010| ) and Beerma n~et al.| ( |2012| ), the sec- 
ond, a more accurate but much slower approach, which builds 
probability maps in the age-mass-extinction space by exploring 
all the nodes of the grid and selects the most probable solution 
(see e.g., |Fouesneau & Lancon|2010| ). 

The method presented here also explores parameter proba- 
bility maps; however, by restricting the analysis to the models 
located in the vicinity of the observed absolute magnitudes (the 
distance to the object has to be known), we significantly reduce 
the computation time. The scheme of the method is sketched in 

A 3D grid of cluster models is built for every value of 
the three physical parameters t, m, E(B - V)Q for simplicity 
Fig.[T](a) shows only a grid for age and mass. For the description 
of the grid, see §|3] Each node of the grid contains 1 000 models 
of the same age, mass, and extinction. They populate the pho- 
tometric parameter space (absolute UBVRI magnitudes). Figure 
[T](b) shows Mu and M B with only 100 models per node without 
extinction, and the used model grid is much more continuous in 
photometric parameter space. 

When the observations are considered in Fig.[T](c), along 
with their error bars (cr; hereafter we use cr = 0.05 mag for all 
the passbands of artificial and real cluster samples studied in this 
paper), which in general can be different for every magnitude, all 
the models situated within 3-cr from the observed magnitudes 
are selected. Fig.[T](d) shows the nodes to which the selected 
models are associated. Other nodes do not play any role in the 
derivation of parameters, therefore the speed of the algorithm 
is increased significantly. Finally, the distributions of age and 
mass, displayed in Figs. [T](e, f)> are derived from the selected 
models, as for the extinction, which is not shown in the figure. 

For the selected models (Fig.[T](c) circle) we apply weights 
as follows: for the models located within l-cr from the observed 
magnitude a weight of 0.68 is assigned, for the ones between 
l-cr and 2-cr - 0.28, and for the ones between 2-cr and 3-cr 
- 0.04. The probability density distributions displayed in Figs. 
[T](e, f) are derived by normalizing the total area of each his- 
togram to 1. The solution is taken as the maximum of these ID 
distributions. We compute confidence intervals (error bars) by 
excluding the first and the last 16% of the area in histograms, 
following the method of "central interval" presented in § 2.5.1 
of |Andrae| ( |20T0] ). 
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Fig. 1. Scheme of the method for deriving physical parameters of 
star cluster. Panel (a) shows the grid of models (only age and mass 
are displayed; same for extinction). These models are represented in 
the photometric space (only M v and M B are shown) displayed in panel 
(b; only 100 models per node without extinction are shown). In panel 
(c) all the models situated within 3-cr (circle) from the observed abso- 
lute magnitudes are selected. Panel (d) shows the nodes associated with 
these models. The ID age and mass distributions of the selected models 
are displayed as normalized histograms in panels (e) and (f), used to 
determine most probable values, indicated by solid vertical lines; confi- 
dence intervals are indicated by dashed lines. 



3. The age-mass-extinction grid of models 

To derive physical parameters of star clusters with the method 
described in §|2] a large age-mass-extinction grid of models is 
computed, by applying the algorithm described by Deveikis et al. 
(2008). The stellar masses are generated randomly sampling the 
IMF (corrected for binaries; Kroupa 2001) and their luminosi- 
ties are derived from stellar isochrone of the selected age and 
metallicity of the cluster model. The process is continued until 
the total mass of generated stars reaches the mass of the cluster 
model. Then, taking the distance to the cluster into account, the 
magnitudes are computed using the Johnson-Cousins UBVRI 
photometric system ( |Maiz A pellaniz 2006]). 

For stellar models, we took the PADOVA isochroneqj from 

ire 



Mar igo et aL] ( |2008| ), with corrections by Girard i et al.| ( |2010| ) 
for the TP-AGB phases. The model grid for a single metallicity 
Z = 0.008 contains the following nodes: ages from log(£/yr) 



We refer to extinction as E(B - V) hereafter. 



PADOVA isochrones from "CMD 2.4": http://stev.oapd.inaf.it/cmd 
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= 6.6 to 10.1 in steps of 0.05, masses from log(m/M ) = 2 to 
5 in steps of 0.05. This gives 71 values of age and 61 values 
of mass, and the grid consists of 1 000 models per node, i.e. ~ 
4 x 10 6 stochastic models. To limit the number of models to 
store in computer's memory, extinction is computed when the 
observed cluster is compared with the grid of models. It ranges 
from E(B - V) = to 1 in steps of 0.02, therefore 51 values for 
the extinction. 

To accelerate computation of integrated magnitudes of 
stochastic star cluster models, we defined a threshold in the 
isochrone, under which the total luminosity of fainter stars is 
computed by integration of the stellar luminosity function along 
the isochrone. Above the threshold, which is defined to select 
20% of the most massive stars, the contribution of the high-mass 
stars is modeled by the algorithm of Devei kis et al.| ( |2Qd8] ). The 
models built by this improved procedure share the same prop- 
erties as models built by simulating all stars of the cluster, but 
require much less computation with a speed gain of a factor -10. 



4. Test of the method on an artificial cluster sample 

4.1. Degeneracies in an artificial cluster sample 

We simulated artificial star cluster samples with known age, 
mass, and extinction and used them as input clusters to evalu- 
ate the ability of our method to derive physical parameters. The 
artificial cluster samples consist of 10 000 clusters with a random 
age in the log(7/yr) range of [6.6, 10.1]. To simulate the mass of 
input clusters, we used a power-law cluster mass function with 
index -2 in the range of log(m/M ) = [2.7, 4.3], so as to have 
more low-mass clusters in the sample. We have two artificial 
samples: one without extinction and the other with E(B - V) in 
the range of [0.0, 1.0] using th e Milky Way standard extinction 
law from Car delli et alJ ( [T9§9] ). 

Figure [2] displays the results for artificial cluster sample built 
without extinction. Panels (a) and (b) show results of the age and 
mass of the artificial cluster sample without introducing pho- 
tometric observation errors. Panels (c) and (d) show the case 
with Gaussian photometric errors of 0.05 mag randomly added 
to each magnitude of the sample clusters. The introduction of 
photometric errors results in broadening around the one-to-one 
line. 

In Fig.|2](b), the asymmetry observed in the mass derivation 
slightly favors high masses. This is because that in the grid of 
star cluster models, the nodes of models with higher mass have 
magnitudes that are less dispersed than the nodes with lower 
mass models, as a consequence of stochasticity (see e.g., |De-| 
|veikis et al.|20 08). Thus, when the models of two nodes of dif- 
ferent masses are located in the UBVRI "sphere" around the ob- 
servation shown in Fig.[T](c), the more massive node, with less 
dispersed models, will dominate. To balance out the effect, a 
cluster mass function could be used to decrease the importance 
of the nodes of massive cluster models. Although Fig.|2](b) 
shows asymmetry, the means and standard deviations given at 
log(m/M ) = 3.0, 3.5, and 4.0 indicate that this phenomenon is 
slight. 

Figure [3] gives the results of artificial clusters in the case of 
extinction. Panels (a), (b), and (c) show the age, mass and ex- 
tinction derived when there is no photometric errors, and pan- 
els (d), (e), and (f) when there are Gaussian photometric errors 
of 0.05 mag included. In Fig.|3](d), in the case of photomet- 
ric errors and extinction, we observe that a broadening around 
the one-to-one relation increases, especially for log(7/yr) > 8, 
which is associated with broadening in the extinction (panel f), 
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Fig. 2. Derived parameters of 10000 artificial clusters without ex- 
tinction. Panels (a) and (b) show age and mass for a sample without 
photometric errors; (c) and (d) for a sample with Gaussian photomet- 
ric errors of 0.05 mag. Panel (b) also shows the means and standard 
deviations (circles with error bars) of the derived mass distributions for 
clusters with true mass log(m/M ) = 3.0, 3.5 and 4.0. The density scale 
is logarithmic. 



which is a sign of the age-extinction degeneracy. It creates two 
streaks above and below the one-to-one relation in the range of 
8 < log(r/yr) < 9.5 that were already perceptible in the case 
without photometric errors in Fig.[3](a). However, including of 
photometric errors does not significantly affect the derivation of 
mass (panels b and e). We note that a gap in derived ages at 
log(£/yr) = 9.15 is a feature of isochrone due to the increase in 
the production rate of AGB stars, which was discussed in Girardi 
|&Bertelli| ( [T998) . 

These degeneracy streaks have first been reported by |Foues 



neau & Lancon ( 2010 ). The degeneracy streak above the one-to- 



one line seen in Fig.[3j(d) concerns clusters that are young and 
that possess intrinsically high extinction, but are derived by the 
method as older and having lower extinction. Conversely, the 
streak of clusters below the one-to-one relation involves objects 
that are derived as older and that have lower extinction than they 
do in reality. These features are important to keep in mind when 
deriving the physical parameters of unresolved star clusters. 

The treatment of the metallicity effects on the derivation of 
age, mass, and extinction of star clusters is the subject of a forth- 
coming paper. To illustrate the effect, we derived the parameters 
of a cluster sample of Z = 0.008 metallicity, successively with a 
model grid of much lower metallicity, Z = 0.00013, and another 
one with much higher metallicity, Z = 0.03. In the first case, 
the ages are systematically overestimated by -0.5 dex, the up- 
per streak is more developed, and the lower one disappears. The 
masses are also overestimated by ~0.5 dex, and there is a pref- 
erence for overestimated extinction. In the second case, the ages 
are only slightly (-0.1 dex) underestimated. The upper streak 
decreases without vanishing, and the lower one becomes more 
populated. The masses are underestimated of -0. 1 dex, and the 
extinction is almost unaffected. 



Article number, page 3 of 



10 




7 8 9 

log(£/yr) true 



log(A7?/M ) true 



0.4 0.8 
E(B-V) true 



Fig. 3. 

(c) show 
mag. 



Derived parameters of 10000 artificial clusters with extinction randomly chosen in the range of E{B - V) = [0, 1]. Panels (a), (b) and 
age, mass, and extinction for a sample without photometric errors; (d), (e), and (f) for a sample with Gaussian photometric errors of 0.05 



4.2. Is it possible to reduce the age-extinction degeneracy? 

The upper and lower streaks in Fig.[3](d) suggest that if a wide 
extinction range is allowed in a simulated sample, then there are 
possibilities that a cluster mimics an older one with lower ex- 
tinction, or inversely a younger one with higher extinction. If 
the true extinction range of the cluster population is narrow, then 
we could restrict the search for the extinction within a narrow 
range in the model grid, resulting in decrease of age-extinction 
degeneracy. 

In Fig. [4] we show results of the tests for a sample of 10000 
artificial clusters with true extinction from a range of E(B - V) 
= [0, 0.5] and with Gaussian photometric errors of 0.05 mag 
randomly added to each magnitude of the sample clusters. This 
cluster sample was studied twice, with different allowed extinc- 
tion ranges in the model grid. 

In the first test, the extinction of the model grid was allowed 
to vary in a wide range, E(B - V) = [0, 1], shown in Fig.|4](a, 
b, c). If a cluster has a true extinction of 0.5, then the maximum 
underestimation of its extinction can be from to 0.5. But if a 
cluster has true extinction of 0, the maximum overestimation of 
its extinction could range from to 1. This explains why the 
lower streak (i.e. clusters with overestimated extinction) is more 
extended than the upper one in panel (a). 

In the second test, the allowed extinction range of the model 
grid was reduced to a range of [0, 0.5]; Fig.|4](d, e, f). The 
constraint on the extinction range resulted in a reduction of the 
lower degeneracy streak, seen in Fig.[4](d), which is less devel- 
oped than in Fig.|4](a). From the comparison of Figs.|4](b) and 
(e), we see that the mass is less affected by the degeneracy. We 
note that only the lower streak is modified, since only the higher 
limit of the allowed extinction range was changed from 1.0 to 



0.5. The upper streak is not modified, because the lower limit of 
the allowed extinction range was not changed. 

To quantify the reduction of the lower degeneracy streak due 
to reduction of the extinction range, all the models situated in the 
square shown in Figs.[4](a) and (d) were counted. There are -480 
clusters in that region for the wide extinction range (panel a) and 
only ~1 10 for the low extinction range (panel d). The reduction 
of the extinction range from E(B - V) = [0, 1] to [0, 0.5] thus 
decreases more than four times the strength of the degeneracy 
streak. 

We conclude that to reduce the degeneracy streaks seen in 
Fig.|3](d), we should make a reasonable assumption on the extent 
of the possible extinction range within the galaxy hosting the 
studied cluster population. 



5. Application of the method to real star clusters 

5.1. The M31 galaxy star cluster sample 

A star cluster catalog of 285 objects located in the south west 
field of the M3 1 galaxy was compiled by Narb utis et al.| ( |2008| ) 
using the deep BVRI and H a photometry from the Subaru tele- 
scope, as well as multiband maps based on HST, Galex, Spitzer, 
and 2MASS imaging. The UBVRI photometry was derived us- 
ing the Local Group Galaxy Survey data from Massey et al. 
( 2006 ). The magnitude limit of the cluster sample was set to 

the V = 20.5 mag. 

Vansevicius et al. ( 2009) selected 238 clusters from that sam- 
ple, excluding the ones with strong H a emission, and compared 
their multiband colors to PEGASE ( |Fioc & Rocca-Vol merange 
1997| ) SSP models to derive their age, mass, metallicity, and ex- 
tinction. Spitzer data was used to constrain the maximum ex- 
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Fig. 4. Derived parameters of 10000 artificial clusters with true extinction randomly chosen in the range of E(B - V) = [0, 0.5] with Gaussian 
photometric errors of 0.05 mag. Panels (a), (b), and (c) show age, mass and extinction when a sample is studied using a grid of models with an 
allowed extinction in the range of [0, 1]; (d), (e), and (f) show results for the same cluster sample studied with a narrower allowed extinction in the 
range of [0, 0.5]. Panels (c) and (f) are denser than in Fig.[3] because the same number of clusters are located on a smaller range of extinction. A 
square indicated in panels (a) and (d) is used to quantify the degeneracy effect. 



tinction for each cluster. They reported ~30 classical globular 
clusters with low metallicity older than 3 Gyr. The remaining 
-210 younger clusters were classified as objects belonging to 
the disk, with average metallicity Z = 0.008. Figure [5] shows the 
U - B vs B - V and U - V vs R- I diagrams of the 238 star 
clusters from Vansevicius et al. (2009), compared to the grid of 
star cluster models built in §[3J 

Since our grid of models has single metallicity, Z = 0.008, 
we attempted to exclude more metal-rich clusters from |Van-| 
sevicius et al.| ([2009 ) sample by only selecting objects with 
galactocentric distance over 7 kpc. This subsample consists of 
216 clusters and is displayed in Fig.[5] It is studied with our 
method described in §[2} usin g the Milky Way standard extinc- 
tion law (Cardelli et al. 1989) and distance modulus to M31 of 



McConnachie et al.| ( [2005] ). 



(m - Af)o = 24.47 derived by 

Figure[6]presents the age (panels a, b, c), mass (panels d, e, f) 
and extinction (g, h, i) of 211 clusters (from the 216 sample) de- 
rived by using our method; for five clusters, no model was found 
within the 3-cr around the observed magnitudes. Clusters were 
studied twice, first with a narrow extinction range E(B - V) = 
[0.04, 0.5] allowed in the model grid, and second with a wider 
one, [0.04, 1.0]. The first column (panels a, d, g) compares 
the age, mass, and extinction derived when a narrow extinction 
range is allowed vs the [Vansevicius et al.| ([2009) results. The 
second column (panels b, e, h) compares the age, mass, and ex- 
tinction derived when a wide extinction range is allowed vs the 
Vansevicius et al. (2009) results. The last column (panels c, f, 
i) compares the results obtained with a wide extinction range 
allowed vs the ones obtained with a narrow extinction range al- 
lowed. 



In Fig.[6](a), when clusters are studied in a narrow allowed 
extinction range, E(B - V) = [0.04, 0.5], the derived ages show 
the same features as the models studied in Fig.[3](d), reproduced 
here in the background of the panel. The degeneracy streaks de- 
velop above and below the one-to-one line, perpendicularly to it, 
and are marked by ellipses numbered "1" and "2" in Fig.[6](a). 
As for the models in the background, the upper degeneracy 
streak concerns clusters with overestimated age and underesti- 
mated extinction, shown by ellipse "1" in Fig. 0(g). In contrast, 
for the lower degeneracy streak (ellipse "2"), the age is underes- 
timated and extinction is overestimated. Since the upper streak 
(ellipse "1") is more developed than the lower one ("2"), we in- 
terpret that as a result of too narrow an extinction range in the 
model grid allowed E(B -V) = [0.04, 0.5]. 

In Fig.[6](b), when clusters are studied in a wider model ex- 
tinction range, E(B-V) = [0.04, 1], the upper degeneracy streak 
is less populated, and the lower one extends, meaning that some 
of the clusters have overestimated extinction (ellipse "3"), also 
shown in panel (h). 

If the intrinsic range of extinction for a cluster sample is 
wide, as in the galaxy M31, then a narrow extinction range of 
models produces extended an upper streak (ellipse "1") and a 
smaller lower streak (ellipse "2"), shown in Fig.[6](a). In panel 
(b), when the allowed extinction range of models is wide, the up- 
per streak retracts and the lower streak develops. We conclude 
that in a galaxy with wide extinction range, it is not possible to 
derive parameters for the clusters affected by age-extinction de- 
generacy and additional constraints for the extinction are needed; 
e.g., [Vanseviciu s et al.| ( |2009| ) used a Spitzer emission map, to 
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Fig. 5. M31 galaxy star cluster sample of 238 objects from Vansevicius et al. (2009 ) (circles and stars) and the studied sample of 216 objects 
(stars) having galactocentric distance over 7 kpc. Panels show: (a) U - B vs B - V and (b) U - V vs R - I diagrams. The model grid (without 
extinction) used to derive physical parameters of clusters is displayed in the background with indicated reddening vectors following the Galactic 
extinction law. The continuous line traces PADOVA SSP of Z = 0.008. 



trace the dust lanes of M3 1 and reduce the age-extinction degen- 
eracy. 

Figs. [6](d, e) show that for high-mass clusters we obtain 



lower masses than given by |Vansevicius et al.|(|2009|). I nspect- 
ing the metallicity values provided by Vanseviciu s et al.| ([2009), 



we note that high-mass clusters have metallicities lower than the 
Z = 0.008 used in our model grid. This could be a sign of the 
metallicity effect on the derivation of physical parameters, and 
will be investigated in the forthcoming paper. 



5.2. The M33 galaxy star cluster sample 



|San Roman et al.| ( |2009| ) reports observing of 161 clusters us- 
ing the HST, allowing them to evaluate their ages and extinc- 
tions by PADOVA [Girardi et al.](|2002| ) isochrone fit to their re- 
solved CMDs. Recently , |Ma|(|2012|) has used 392 clus ters from 
the compiled catalog of Sarajedini & Mancone (2007) and im- 
ages from |Massey et al.| ( |2006| ) to provide UBVRI broad-band 
integrated photometry data. 

There are 40 clusters common to both the San Roman et al. 
( 2009 ) and the |Ma| ( |2012[ ) catalogs, making them interesting 
to compare the paramete rs derived from the res olved method 
(isochrone fit to CMDs by San Roman et al. ( 2009 )) and stochas- 
tic method, which was applied here to the integrated photometry 
data of [Ma1 ( |20T2l ). 

Figure |7](a) presents the U - B vs B - V diagram of the 392 
clusters from the Ma (2012 ) catalog compared to the grid of clus- 
ter models. The selected 40 star clusters are also in dicated. To 
account for calibration of the B band, as discussed by Ma (2 012|), 
where he compared his photometry to previous wo rks by|Sara- 



jedini & Mancone| ( [2007] ), |Park & Lee| ( [2007] ), and |San Roman 



et al. ( 2009| ), we shifted the B band to brighter magnitudes by 
0.1 mag for all clusters. We note that the parameter derivation 
results do not change when the B band is not taken into account. 



Since the metallicity content of the M33 galaxy is similar to 
the one of the LMC ( [Bresolin et al.||2010|), we assume tha t the 
interstellar extinction law derived by |Gordon et al.| ( |2003| ) for 
the LMC can be applied to the M33 cluster sample. As in |San 
Rom an et al.| ( |2009|), we adopted the M33 di stance modulus of 



(m - Af)o = 24.54 ( [McConnachie et aTp'0 05 ). 

For the M33, |U et al.| ( [20091 ) derived the extinction of 22 
supergiants, which resulted in an extinction distribution cen- 
tered on E(B - V) = 0.1. |U et al.| ([20091) also used the data 



of [Rosolows ky & Simon| ( [2008 ) to derive E(B - V) values for 
58 HII regions, and show in their Fig. 9 that extinction can be 
expected to be E(B - V) < 0.3 for those regions (except for 3 ob- 
jects), with an average E(B - V) ~ 0.11. A foreground Galactic 
line-of- sight extinction in the direction of M33 of E(B - V) = 
0.04 is estimated from Schlegel et al.| ( |1998] ) extinction maps. 

We studied the M33 clusters using two different allowed ex- 
tinction ranges for the model grid. The first one is narrow, ex- 
tending between the foreground Galactic line-of- sight extin ction 
in the direction of M33 up to the higher limit given by |U et aL 
( [2009] ), i.e, E(B - V) = [0.04, 0.30]. The second one is a wider 
extinction range, E(B-V) = [0.04, 1.0], such as the one used for 
the M31 case. 

Figure [8] presents the age (panels a, b, c), mass (panels d, 
e, f), and extinction (g, h, i) of the 40 clusters derived by using 
our method. The results obtained in the narrow extinction range 
(panels a, d, g) and the wide extinction range (panels b, e, h) are 
compared to the ones obtained by San Rom an et"aL] ( |2009| ) using 
isochrone fit to CMDs. Panels (c, f, i) compare the results we 
obtained with the wide extinction range allowed vs the narrow 
extinction range allowed. It shows that the difference is relatively 
small between the two kinds of solution, showing that the age- 
degeneracy is playing a minor role in this M33 cluster sample. 

The derived mass does not differ significantly, as seen in 
Fig.[8](f). Panels (g, h, i) reveal that only one fourth of the clus- 
ters are affected by an increase in extinction when the allowed 
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extinction range is widened. Only some of them have a much 
higher derived extinction when the wide allowed extinction 
range is used. These few clusters affected by the age-extinction 
degeneracy can be seen in panel (c), appearing younger when a 
wide extinction range is allowed than when the range is narrow, 
as already seen in M31 cluster sample in Fig.[6](c). 

6. Conclusions 

We presented a method that aims to derive the age, mass, and ex- 
tinction of unresolved star clusters by using broad-band photom- 
etry and taking the stochastic sampling of stellar masses in clus- 
ters into account. We investigated the behavior of the method 
on a sample of artificial star clusters, in order to trace the dif- 
ferent degeneracies between parameters, for different choices of 
photometric errors and extinction. 

These tests allowed us to quantify the age-extinction degen- 
eracy. We demonstrated that for the intrinsically narrow extinc- 
tion range of the star cluster sample, the age-extinction degener- 
acy can be resolved even in cases where individual exact extinc- 
tion values are not known for each cluster. 

The age-extinction degeneracy has been observed in the real 
star cluster sample of Vanseviciu s et al.| ( |2009| ) for M3 1 and to 
a lesser ex tent in the one c omposed of clusters common to the 
|Ma| ( [2QT2| ) and |San Roman et al.| ( [2QQ9l ) catalogs for M33. The 
physical parameters derived by our method for different extinc- 
tion ranges in each case, have been compared to the values pro- 
vided in these studies, showing the impact of the age-extinction 
degeneracy, especially when the true extinction range of the clus- 
ter population is wide. 

The M31 star cluster sample from|Yansevici us~et al.| (f2QQ9 ) 
showed that a true extinction range in this galaxy is wide enough, 
so that the age-extinction degeneracy is significantly developed, 
making the parameters difficult to derive, especially for older 
ages. In such cases, it is preferable to use external constraints on 
the extinction of individual clusters. 

The M33 star cluster sample shows little sign of age- 
extinction degeneracy, since there is only a small difference be- 
tween solutions in narrow and wide extinction ranges. The range 
of allowed extinction of E(B - V) = [0.04, 0.30] gives parame- 
ters consistent with the results of isochrone fit to CMDs by [San] 
|Roman et al.| ( |2009] ). 

The follow-up paper will be dedicated to study of the metal- 
licity effects to derive physical parameters of star clusters and 
to gain a complete understanding of degeneracies introduced by 
stochastic effects. 
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Fig. 6. M31: age, mass, and extinction of 211 star clusters from ( Vansevicius et al.|2009 as V09). Panels (a, d, g) show results when a narrow 
extinction range of E(B - V) = [0.04, 0.5] is allowed in the model grid, compared to V09 values. Panels (b, e, h) show results for a wider allowed 
extinction range [0.04, 1.0], compared to V09 values. Panels (c, f, i) compare results obtained in a wide extinction range vs the ones obtained in a 
narrow extinction range. In panels (a, b, c), a density map of Fig.[3](d) is reproduced to show regions of age-extinction degeneracy. Error bars are 
computed as described in § 2. Six clusters from V09, which have E(B - V) > 1, are not shown in panels (g) and (h). Dashed ellipses, numbered 
"1", "2", and "3", describe majority of associated clusters in the age and extinction panels. 
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Fig. 7. M33 galaxy star cluster sam ple of 161 objects from cata log of San Roman et al.|f2010| ) (circles and stars) and the studied sample of 40 
objects (stars) common to Ma ( 2012 ) and San Roman et al. ('2009). Panels show (a) U - B vs B - V, and (b) U - V vs R - 1 diagrams. The model 
grid (without extinction) used to derive physical parameters of clusters is displayed in the background with indicated reddening vectors following 
the LMC extinction law. The continuous line traces PADOVA SSP of Z = 0.008. The B band magnitudes of clusters were shifted by 0.1 mag. 
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Fig. 8. M33: age, mass and extinction of 40 star clusters common to catalogs of Ma ( 2012| ) and ( |San Roman et al. 2009 as S09). Panels (a, d, g) 
show results when a narrow extinction range of E(B - V) = [0.04, 0.3] is allowed in the model grid, compared to S09 values. Panels (b, e, h) show 
results for a wider allowed extinction range of E(B - V) = [0.04, 1.0], compared to S09 values. Panels (c, f, i) compare results obtained in a wide 
extinction range vs the ones obtained in a narrow extinction range. In panels (a, b, c), a density map of Fig.[3](d) is reproduced to show regions of 
age-extinction degeneracy. Error-bars are computed as described in § 2. 



Article number, page 10 of 10 



